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Abstract 

The high-barrier quantum tunnehng regime of a Bose-Einstein condensate confined in a ring- 
shaped optical lattice is investigated. By means of a change of basis transformation, connecting 



'^ ■ the set of 'vortex' Bloch states and a Wannier-like set of localized wave functions, we derive 

O . a generalized Bose-Hubbard Hamiltonian. In addition to the usual hopping rate terms, such a 



Hamiltonian takes into account interaction-driven tunneling processes, which are shown to play a 



^ ' principal role at high filling factors, when the standard hopping rate parameter turns out to be 

00 ■ 
00 . negative. By calculating the energy and atomic current of a Bloch state, we show that such a 

hopping rate must be replaced by an effective hopping rate parameter containing the additional 
^^ '_ contribution an interaction-driven hopping rate. Such a contribution turns out to be crucial at 

high filling factors, since it preserves the positivity of the effective hopping rate parameter. Level 
- ^-. , crossings between the energies per particle of a Wannier-like state and the superfluid ground state 

X, 

H are interpreted as a signature of the transition to configurations with macroscopically occupied 

states at each lattice site. 

PACS numbers: 03.75.Lm, 03.75.Hh, 03.75.Kk 



I. INTRODUCTION 

The study in the last decade of ultra-cold bosonic atoms in optical lattices has enabled 
the realization of an active and fruitful convergence of atomic and condensed matter physics. 
Particularly, the analogy of such systems with a solid material, where the bosons play the role 
of the superconducting electron pairs and the laser beams act as the ionic crystal, became 
the leitmotiv of numerous applications (l|, |2|. In their seminal experiment, Greiner et al. [3|] 
showed that by increasing the lattice potential depth in a three-dimensional optical lattice, a 
quantum phase transition from a superfluid state to a Mott insulating state can be achieved. 
This had been predicted by Jaksch et al. [4], who accurately described such a transition 
within a Bose-Hubbard model at filling factors of the order of unity. Actually, most research 
has so far been focused on optical lattices with such a low filling factor, whereas the high 
filling factor domain appears scarcely treated. Such high-filling configurations are expected 
to be noticeably affected by the on-site interaction between bosons, as the Wannier single- 
particle ground-state wave function in every site should be replaced by a macroscopic wave 
function [5] . A suitable configuration to experimentally investigate this type of condensates 
could be given by a ring-shaped lattice, where a toroidal trap becomes symmetrically divided 
by a number of potential barriers radiating away from the trap center [6|]. In fact, apart 
from presenting the ideal geometry to sustain persistent currents, such a lattice would also 
exhibit a perfect azimuthal periodicity for any number of lattice sites. This would permit 
to achieve extremely high filling factors within the present experimental possibilities for the 
maximum number of particles in the whole condensate. The effect of raising a single barrier 
across a long-lived persistent current in a toroidal condensate, has recently been investigated 
as the first realization of an elementary closed-loop atom circuit [7]. The generalization of 
such experiments to ring lattices has shown to be quite attainable in the light of the works 

n n 

of Amico et al. [8] and Henderson et al. [9]. In fact, while the former have thoroughly 
discussed the experimental setup for realizing a ring lattice, such a system was actually 
generated by the latter, utilizing a rapidly moving laser beam that 'paints' a time-averaged 
optical dipole potential, transforming a toroidal condensate into a ring lattice. 

From a theoretical viewpoint, recent investigations have analyzed the effect of rotation 
on the ground state properties of bosonic atoms confined in a one-dimensional ring lattice 



at low filling factors 10|]. A nonrotating ring lattice, on the other hand, has been predicted 



to sustain persistent currents [ll| if the phase difference between adjacent sites takes certain 



values 



U 



In addition, the buildup of winding number in the phase transition from Mott 
insulator to superfluid driven by tunneling rate increase, has been shown to proceed through 
the so-called Kibble-Zurek mechanism, except for very slow quench times jl2| . In the present 
work we will concentrate our attention on such nonrotating configurations with high barriers 
and high filling factors. The starting point of a theoretical approach to this kind of systems 
should consist in exploring an adequate variant of the Bose-Hubbard (BH) model, which 
should be expected to exhibit occupation dependent parameters [l3|, llj]. As usual, the main 
ingredient to derive such a BH Hamiltonian consists in finding a suitable set of orthogonal 
Wannier-like functions, for which a number of variational schemes have been proposed |14J ~ 



18|. Here, rather than resorting to such methods, we shall obtain our set of Wannier-like 
functions simp ly as a 'basis change' from the orthogonal set of stationary 'vortex' Bloch 
states 



19 



20|]. Then, it will be shown that such functions possess the main properties of 
the single-particle Wannier functions, except for their dependence on the filling factor, and 
thus they become the adequate tool to study the slightly perturbed Bloch states arising 
from small occupation number imbalances, or from small changes on the relative phase 
between adjacent sites. Under such conditions, a generalized BH Hamiltonian that takes 
into account interaction-driven tunneling processes will be derived. Such contributions. 



which were previously investigated for double- and triple- well configurations 2l|-|23j, will 
be shown to play a principal role at high filling factors. Finally, by considering the level 
crossing between the energies per particle of a Wannier-like state and the superfluid ground 
state, we will discuss the transition to configurations with a macroscopic occupation at each 

site. 



This paper is organized as follows. In Sec. II, we describe the ring lattice and remaining 
condensate parameters. In Sec. HI, we analyze the main properties of Bloch and Wannier- 
like states. In Sec. IV, we derive the generalized BH Hamiltonian from which the energies 
of Bloch states are calculated, and the continuity equation at a given lattice site and the 
corresponding atomic current are extracted. Finally, in Section V we discuss our numerical 
results for the tunneling parameters and level crossings, while in Sec. VI we present our 
summary and main conclusions. 



II. RING-SHAPED LATTICE AND CONDENSATE PARAMETERS 

We consider a Bose- Einstein condensate of rubidium atoms confined by an external trap 
l^rap, consisting of a superposition of a toroidal term VJ;oro and a lattice potential Vi^ formed 



2J,|25|, the toroidal 



by radial barriers. Similarly to the trap utilized in recent experiments 
trapping potential in cylindrical coordinates reads, 

Koro(r, ^) = Y [ooy + ooy] + Vo exp(-2 rV A^) (1) 

where Ur and Uz denote the radial and axial frequencies, respectively, and M denotes the 
atom mass. We have set u^ » uJr to suppress excitation in the z direction. In particular, we 
have chosen u;j./(27r) = 7.8 Hz and ujz/{2'n-) = 173 Hz, while for the laser beam we have set 



Vq = lOOhur and Aq = 6/,-, with Ij. = Jh/{MiOr)- On the other hand, the lattice potential 
is formed by N^ Gaussian barriers of width A^ and amplitude Vb, located at equally spaced 
angular positions Ok = ^Tik/N^, where -[[{N^ - l)/2]] < A; < [[A^c/2]] with [[■]] denoting the 
integer part, 

Vi^[x, y) = Vb 2^ 6[sm(6'fc) y + cosiOk) x\ exp <^ ^ > , (2) 

where denotes the Heaviside function. 

In the mean-field approximation, the stationary states are solutions of the Gross- 



Pitaevskii (GP) equation [26| 



-^V^ + Krap(r) + ^iV|^(r)|2| V'(r) = /xV^(r) (3) 

where N, fi and V'(r) respectively denote the number of particles, the chemical potential and 
a two-dimensional (2D) order parameter normalized to one [27|. The effective 2D coupling 



constant g = g-^£)J Muz/^nfi is written in terms of the 3D coupling constant between the 
atoms g^o = 4:7cah /M, where a = 98.98 oq denotes the s-wave scattering length of ^'^Rb, Uq 
being the Bohr radius. 

III. BLOCK AND WANNIER-LIKE STATES 

We shall restrict our treatment to the case of high enough barrier heights, where quantum 
tunneling between sites turns out to be the dominant dynamical process. Such a regime 



arises when the ground-state chemical potential becomes smaller than the minimum of the 
effective potential barrier dividing two lattice sites IgI . The most general solution of the GP 
equation (j3]) is given by a Bloch state of the form [l9|, |20| 



ipm{r,9) 



Amd 



fm{r,0), 



(4) 



where /„ 



is invariant under rotations in 211 /N^ and the winding number m plays the role 



of an 'angular' pseudomomentum satisfying the constraint — [[(A^c — l)/2]] < m < [[Nc/2]]. 
Such a constraint arises from the fact that all possible solutions can be reduced to those 
existing in the first Brillouin zone in pseudomomentum space 20|]. We shall restrict ourselves 
to Bloch states of the lowest energy, i.e., to the ground 'vortex' states [6]. In the language of 
crystal lattices, we would say that we shall restrict our treatment to the subspace of Bloch 
states of the 'ground band'. In such a context, the orthogonality of a pair of Bloch states, 
ipm and tpn, can be easily proven as follows. First, the corresponding integral may be split 
into separate integrals over each site, where we make the change of variable 6' = 9 — 9^. 
Then, taking into account the rotational symmetry of the corresponding functions fm and 
/„, along with the equality 



^exp[i(m - n)9k] = 5m,nNc 

k 



(5) 



the orthogonality can be demonstrated. 

Now, taking into account the periodicity of a Bloch state in the reciprocal lattice, 
ipm+jNc = i^m, it must havc a Fourier series expansion with 'wave vectors' 9'j. = {9k + 
6'fc_|_i)/2 = 9k + tt/Nc in the direct lattice as follows 28| 29|, 

1 



i^m{r,9) 



^ 



J2'^kir,i 



^idf^m 



(6) 



c k 



where the Fourier coefficients in ([6]) are given by the inversion formula 30 1 

1 



Wk{r,9) 



^r 



5:^„(r,^)e-^^ 



(7) 



with the summation over the angular pseudomomentum n being restricted to the first Bril- 
louin zone. Replacing (Hj) in ((Tj) and taking into account the symmetry of /„, we may realize 
that the Fourier coefficients arise from a single function w(r, 9) as follows, 

1 



Wk{r,9) =w{r,9 -9k) 



^ 



Y.fnir,9-9k)e 



in{e-e^) 



(8) 



Thus, pushing forward with the analogy to crystal lattices, we could name the function 

w{r,e) = w,{r,9) = ^^^„(r,0), (9) 



the 'Wannier' function of the ground band ! 



'c n 



Although we shall see that it shares many 



formal properties with the well-known Wannier functions, we shall also show that it presents 
a remarkable difference. So, we feel more appropriate to speak in the following of a Wannier- 
like function. Let us first show the similarities. Taking into account that the Bloch 'vortex' 



states |20] fulfill ip^ = ip-n, while the ground and highest states, ipo and ipNc/2 {Nc even), 
respectively, are real, it is easy to show that w{r,6) given by ([9]) must be a real function. 
Also, from the orthonormality of Bloch wave functions and Eq. ([5]), one may readily check 
that the set of Wannier-like functions centered on different fc-sites, w{r, 9 — 6*^), form indeed 
an orthonormal basis of the subspace of Bloch states of the ground band. In addition, given 
that the Bloch 'vortex' states fulfill V^nl^, —&) = V'nl^) ^)j i^ is easy to show that w{r, 9) turns 
out to be an even function of 9 for odd N^.. On the other hand, by considering the rotation 
of the coordinate system in tt/Nc shown in Fig. [H which makes the Bloch wave function 
ipNc/2{r, 6*) an even function of 9, the same parity property may be readily extended to the 
case of Nc even. Finally, by replacing the numerical solutions of the GP equation, ipn{f, &)-, in 
Eq. ([9]), we have shown that our Wannier-like function is indeed a well-localized one, as seen 
in Fig. [H However, there is a most remarkable difference between such a localized function 
and a 'true' Wannier function, which consists in that only the former turns out to depend 
on the filling factor, i.e. the average number of particles at each site, as clearly observed in 
Fig. [H In fact, only for noninteracting bosons our Wannier-like function would not depend 
on the filling factor. We have performed a calculation of the overlap between the Wannier- 
like function given by ([9]) and the corresponding Wannier-like function for noninteracting 
bosons, i.e. with a vanishing coupling constant (? = 0, which yielded 1.00, 0.98 and 0.48, 
for filling factors 5, 62.5 and 6250, respectively. Particularly, the last two values correspond 
to the filling factors of the top and bottom panels of Fig. 1, respectively. Therefore, we 
may conclude that only for filling factors below ~ 60, the Wannier-like functions should be 
almost independent of the average occupation number. We will have more to say about this 
filling factor dependence in the following Sections. 
To conclude it is instructive to rewrite Eq. ([6]) as 

^^(r,^) = -^^^(r,^-e,,)e^™'N (10) 
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FIG. 1. Isocontours of the ground-state wave function density iV'oP (left panels) and of the Wannier- 
hke function density w'^ (right panels). The rotation of the coordinate system denoted by dashed 
lines in the top-right panel makes the Wannier-like function symmetric with respect to the angular 
variable 9. The condensate parameters are Vb/hur = 10, N = 10'^ (top), and V}j/hujr = 80, A^ = 10^ 
(bottom), while the number of lattice sites and the Gaussian barrier width are given by N^ = 16 
and Xh/lr = 0.5, respectively. 



and notice that the above representation will be accurate to the extent that each site presents 
an almost uniform phase, which is consistent with a tight-binding scenario of high barriers 
with a low particle current |6|. 



IV. BOSE-HUBBARD MODEL 



The above similarities between the Wannier-hke function and the 'true' Wannier func- 
tions, offer the adequate framework to estabhsh a BH model for our ring-shaped optical 
lattice. As usual |2|, y, |3l|, the starting point is the second-quantized Hamiltonian 



H 



d^riiHr) 



2M 



V' + Vt 



trap^ 



^(r) + - fdh ^"^(r)^"^(r)^(r)^(r), (11) 



where \l/(r) is the boson field operator. We are interested in slightly perturbed Bloch states, 
which could be given by, e.g., a small relative imbalance between the average population 
of two neighboring sites. Then, for low enough temperatures, such configurations will be 
conveniently described by expanding the field operators in our Wannier-like basis of the 
ground band 

^{r) = ^w{r,e~dk)dk, (12) 

k 

where the operator d^ destroys a particle in the fc- Wannier state and satisfies the usual Rose 
commutation relations. Here we remark that a possible dependence of the operators a^ on 
the filling factor should, under the above conditions, be negligible. In fact, in the previous 
Section we have seen that this is actually the case for filling factors below ~ 60, while for 
higher fillings, only configurations that present small population imbalances should be taken 
into consideration. 

Then, replacing the field operators in flTTl) through Eq. flT2l) and assuming the tight- 
binding limit, where only the coupling to the nearest neighboring states of any given 
Wannier-like state is taken into account, we obtain the following BH Hamiltonian 

HbH = ^ Z^ ^k^k — JZ^i^-k^k+l + dfc+iQ-fc) 



"o" Z^ [o^k^k^kif^k+l + O-fe-l) + (%+l + 0'k-l)^k^kdk 



u 



with 



ak^k^k^k, 



e = d r voir, 



h' 



2M 



V + Krap(r 



J 



d r w(r. 



h' 



2M 



V + Krap(r 



w{r, 6) 
w{r,9±2n/Nc 



(13) 

(14) 
(15) 



J' = -2g (fr w'^ir, 9) w(r, 9 ± 271 /N^) (16) 

U = g I dh w\r,9), (17) 

where the equivalence between the 'ib' expressions at the right-hand side of ( ITSl) stems 
from the reahty of the Wannier-hke functions, while the corresponding equivalence in ( 1T6|l 
results from the parity property of such functions. In addition to the usual tunneling terms 
proportional to the standard hopping rate J, we have also retained in ( 1T3|) interaction 
terms up to the first order in the product of adjacent Wannier-like functions, which are 
proportional to the tunneling parameter J' 32|]. Later we will show that such interaction 
terms may constitute the most significant contribution to the tunneling rate at high filling 
factors. The case of two sites A'"c = 2 is somewhat special since it is the only configuration 
presenting a single neighbor for each site. Then, the expression (TT3ll reduces to 



H, 



BH 



sN 



J 



{N 



N. 



-J' 



(agfli + ajao) + — (aoO-oa^oa-o + alalaiai; 



(18) 



it; 



where, for a fixed number of bosons A^, the particle number operator A^ = agSo + fi-i^i niay 
be replaced by a c-number. Thus, we may see that the only difference with the standard 
two-mode BH Hamiltonian consists in that the standard hopping rate J is replaced by an 
effective hopping rate. 



J. 



eflf 



J 






J', 



(19) 



which includes the additional contribution of an interaction-driven hopping rate ^^J' stem- 



Nc 



ming from boson interactions. Here it is worth noticing that an extended two-mode ap- 
proach, which includes terms in the BH Hamiltonian beyond the present approximation, 



has been recently investigated 21 



33|. 



Next we obtain the mean value of the BH Hamiltonian {N,m\HBH\N,m), where \N,m) 
represents the quantum state of A^ bosons condensed in the Bloch state dl]) of winding 



number m 



19 



20|. To calculate such a matrix element, we may change to the Bloch basis 



in ( 1T3l) by means of the expansion (cf. (^^] 



^r 



E"n 



-indu 



(20) 



where the operator dj, creates a particle in the corresponding Bloch state. Then, a straight- 
forward calculation yields 

{N,m\HBH\N,m)/N = e + ^^^ ^^ | - z/J^ff cos(27rm/Ar,), (21) 



Er, 



where v denotes the number of neighbors [y = 2 [y = V) for Nc > 2 {N^ = 2)). We note 
that the above expression coincides with our previous result [6| in the hmit A^ ^ 1. 

It is instructive to analyze the continuity equation for the /c-th site of a lattice with 

iV, >2, 

d i - 

— (d^Sfc) = — [-f^B/f, a^flfc] = Jk-i^k — Jk^k+i, (22) 

where Jk^k+i denotes the current operator for atoms that move from site k to site k + 1, 

Jk^k+1 — T[(^k+l^cS (^k - O-fc-'eflf «fc+lj; l^-^j 

which has been written in terms of the hopping operator between sites k and k + 1 defined 

by 

T 

Jcs = J + y (4«fc + 4+iafe+i)- (24) 

The mean value of the current operator (1231) for a condensate of A^ particles in the Bloch 
state of winding number m reads 

{N,m\Jk^k+i\N,m) = 2Jeff— sin(27rm/Ar,), (25) 

which does not depend on the site we are considering, as expected. Note that analogously 
to the mean value of the angular momentum [6], the current turns out to be a sinusoidal 
function of the winding number. Note also its proportionality to the effective hopping rate 
Jeff, whereas for the standard BH model such a current turns out to be proportional to the 



standard hopping rate J 



34|. 



The value of the BH model parameters (1141) - (1171) can be easily extracted from the mean- 
field energy of Bloch states Em (see Appendix A). Particularly, from the single value of 
energies of the ground state £q and the highest excited state £n^/2 {Nc even) one obtains 

e = l{£k/2 + So), (26) 

N 

f^ = Tr^(^r/2 + 0, (27) 



1 

2^' 



J = -A/2-£o) (28) 



/ = ^^;(frT)(^;^%-^r) (29) 

where, according to Appendix A, the superscripts 'mt' and '0' denote interacting and non- 
interacting contributions to the energy, respectively. 

10 



According to the Hainiltonian (fT3|) . the energy per particle in a Wannier-hke state, i.e., 
neglecting tunneling processes, is given by 

It is interesting to compare the above energy to the energy per particle of the superfluid 
ground state Eq. Then, from Eqs. f l2T|) and fl30|) we obtain 

Ew-Eo = Z/Jeff - (^^^) |- (31) 

The existence of a superfluid to Mott insulator transition requires the above difference to 
be positive for low barrier heights (superfluid regime), and negative for high barrier heights 
(Mott insulator state). Thus, the level crossing at an intermediate barrier height arising 
from (jH]), should be representing a transition to configurations where the system is well 
described by Nc macroscopically occupied states. We may utilize the above expression to 
obtain the value at such a level crossing, rjcr, of the dimensionless scaling parameter J2|, |3l| 

V = -^, (32) 



relevant to the superfluid to Mott-insulator transition. Thus, assuming Ew = Eq in (13111 . 
we obtain 

r],r = 2Nj{N,-l). (33) 

We may compare the above result with theoretical estimates focusing on critical values of 
the parameter U / J. In fact, it has been pointed out that the superfluid to Mott-insulator 
transition in a one-dimensional BH model can be described by the (1+1)D 0(2) model, 
which gives [l| 

{U/J)cr = 2.2n (34) 

for filling factors n ^ 1 . The above proportionality to the filling factor is also predicted from 
the mean-field Gutzwiller ansatz, which yields {U/J)cr = 2 (v^+ V^ + 1)^ — 8 ra for n ^ 1 
[10|. However, we must recall that mean- field theories only provide a qualitative analysis 
in ID systems. We must also remark that the result (l3l|) arises from a BH Hamiltonian 
that does not take into account the contribution of the interaction-driven tunneling terms 
proportional to J'. So, such an estimate should only be reliable for J^s — J, i.e., for J ^ nJ' 
(?7, :^ 1). However, we shall see in the following Section that such conditions are difficult to 

11 



reach within our condensate parameters. Moreover, we shall show that the standard hopping 
rate J becomes negative above certain filling factor, which means that the parameter U / J 
should actually increase with the average occupation number until becoming divergent and 
meaningless above such a filling factor. 

V. NUMERICAL RESULTS 

In the BH model for linear lattices it is common to measure energies in units of the recoil 
energy E^ = h^k\/2M^ where the Bragg momentum ks corresponds to a lattice potential 
of the form ~ svo?{kBx). To adapt this definition to the present case, we first note that a 
lattice potential ~ sin^(A^c 6'/2) would have the required angular periodicity of l-a jNf,- Then, 
recalling that without barriers the excitation energy per particle of a Bloch state of angular 
pseudomomentum m reads Km"^, where p] 



K = — J -[Mr)Ydr, (35) 

we may realize that our 'recoil energy' should be written 

Er = K{Nj2f. (36) 

We have performed numerical simulations for three particle numbers, A^ = 80, 10^ and 10^; 
given that the corresponding recoil energies turned out to be 0.743 TiUr, 0.740 TiUr and 0.713 
hujri respectively, showing figures that approximate the harmonic energy quantum hujr-, we 
decided, for the sake of simplicity, to keep such a value as our energy unit in all cases. 

We have numerically evaluated the BH parameters through Eqs. fl26j) - fl29j) for the above 
particle numbers and three numbers of lattice sites, N^ = 16, 8 and 4. In Table I, we display 
our numerical estimates for the level crossings for the different condensates and a Gaussian 
barrier width Xb/lr = 0.5. Apart from the dependence of the barrier height parameter VJ,, it 
is interesting to compare the minimum of the effective potential barrier dividing two lattice 
sites Vmin m, with the ground-state chemical potential fi. Recall that in Ref. [6j] we have 
identified the lower bound of the quantum tunneling regime as V^ij^/fi ~ 1. To scale out the 
dependence of the barrier height for different particle numbers, we have represented Vmin and 
/i in units of the chemical potential at zero barrier /iq for each particle number. Compare 
also the numerical estimates for rjcr shown in Table I to those given by the expression (133|) . 

12 



TABLE I. Level crossings arising from Eq. ()3ip . see text for explanation. 



Nc 


N 


N/N, 


Vb/hUr 


Knm//WO 


^/^o 


Vcr 


Vb/huJr \1] 


16 


80 


5 


10.4 


1.29 


1.05 


2.17 


14.1 




103 


62.5 


15.4 


1.39 


1.08 


2.18 


30 




10^ 


6250 


95.4 


1.72 


1.25 


2.23 




8 


80 


10 


4.65 


1.11 


1.01 


2.29 


8.2 




10^ 


125 


10.1 


1.23 


1.03 


2.30 


15.2 




10^ 


12500 


81.4 


1.53 


1.11 


2.32 




4 


80 


20 


2.95 


1.06 


1.00 


2.23 


5.8 




103 


250 


8.93 


1.20 


1.01 


2.68 






10^ 


25000 


77.1 


1.47 


1.05 


2.48 





namely rjcr = 2.13, 2.29 and 2.67 for N, = 16, 8 and 4, respectively. Here it is worthwhile 
pointing out that a similar agreement was found for wider Gaussian barriers {\b/lr = 1). 
Finally, the last column of Table I shows the critical estimate for the barrier height parameter 
Vb arising from Eq. ( IMl) . The absence of data for filling factors above 125 corresponds to the 
negative values obtained for the hopping rate J . Note also that such critical barrier heights 
turn out to be always higher than those of the fourth column, as expected. 

In Figs. [2] and [3] we depict the standard hopping rate J flTSj) computed from Eq. fl28|) . the 
tunneling parameter ^^J' ( TT6|) computed from Eq. (!29|) . and their sum Jefj, as functions of 
the barrier height. Particularly, Fig. |2] shows that while the interaction component ^^J' 



N, 



turns out to be almost negligible for 80 particles and N, = 16 (filling factor = 5), for fewer 
lattice sites it shows a relative increase until becoming of the same order of the hopping rate 
J for Nc = 4: (filling factor = 20). The larger filling factors of A^ = 10^ yield an interaction 
component that turns out to be always larger than the standard hopping rate, to such an 
extent that J eventually becomes negative for N, = 4. Note that such a dramatic change of 
sign occurs in between fillings of 125 and 250 particles (Table I). Negative values of J were 
also predicted by Ananikian et al. for large atom numbers in a double- well condensate 21 1. 
Finally, the extremely high fillings of A^ = 10^ yield again a negative J, as expected, while 
a sort of saturation in the relative weights of the interaction component and the negative 
hopping rate is observed. This is reflected through the quite similar plots of Fig. [3l despite 
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FIG. 2. Standard hopping rate J, interaction-driven hopping rate -^—J', and effective hopping 
rate Jeff, as functions of the barrier height Vf, for the condensates of 80 particles (left) and 10^ 
particles (right). All quantities are given in units of TujO^. In panel (c) the standard hopping rate 
J turns out to be negative for 10^ particles (right), so we have depicted its absolute value \J\. 
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FIG. 3. Absolute value of the standard hopping rate \J\, interaction-driven hopping rate -j^J', 
and effective hopping rate Jeff, as functions of the barrier height V5 for the condensate of 10^ 
particles. All quantities are given in units of fiUr- 
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FIG. 4. Diniensionless scaling parameter rj = (//(vJcs) as a function of the barrier height Vf, for 
Nc = 16 and the condensates of 80, 10^ and 10^ particles. The vertical lines correspond to values 
in the fourth column of Table I. 

the vertical shift for varying Nc, which arises from a decrease of the probabiUty of tunnehng 
events as the number of barriers is lowered. 

To conclude this Section, we display in Fig. H] the diniensionless scaling parameter rj = 
U/luJcs) versus the barrier height, for each number of particles and Nc = 16. Notice that 
quite similar ranges of t] are obtained irrespective of the barrier height interval, and this 
behavior repeats for the remaining values of Nc- 

VI. SUMMARY AND CONCLUDING REMARKS 



We have analyzed the high-barrier quantum tunneling regime of a Bose- Einstein conden- 
sate confined in a ring-shaped optical lattice. Representing the orthogonal set of 'vortex' 
Bloch states through a basis of well-localized Wannier-like functions, we were able to formu- 
late a variant of the Bose-Hubbard model, adequate for slightly perturbed Bloch states at 
any filling factor. In addition to the usual hopping rate terms, such a Hamiltonian contains 
interaction-driven tunneling terms, which are shown to play its most important role when 
the standard hopping rate parameter becomes negative at high filling factors. In fact, by 
calculating the energy and atomic current of a Bloch state, we have shown that the stan- 
dard hopping rate parameter must be replaced by an effective hopping rate containing the 
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additional contribution from the interaction-driven tunneling terms in the BH Hamiltonian. 
We remark the importance of such an interaction-driven hopping rate parameter, since it is 
shown to preserve the positivity of the effective hopping rate at high filling factors. A quite 



similar behavior for such 



two-well configuration 21 1. 



lopping rates was recently predicted for large atom numbers in a 



We have found that, as the barrier height is increased, the energies per particle of a 
Wannier-like state and the condensate ground state exhibit a level crossing, which is in- 
terpreted as a signature of the transition to configurations with macroscopically occupied 
states at each lattice site. It is also shown that the dimensionless scaling parameter, relevant 
to the superfluid to Mott insulator transition, takes a remarkably simple expression at the 
level crossing, which only depends on the number of lattice sites. 

Finally, we would like to point out that a future direction of the present studies will 
consist in exploring the Boson Josephson-j unction dynamics described by a generalized N^- 



mode GP equation |2ll. l22l. l35l|. 
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Appendix A: Alternative calculation of the BH model parameters 

An alternative calculation of the BH model parameters ( TT^ - (TT7j) . which avoids to change 
to a smaller numerical grid to deal with the tiny regions where the integrands of the tunneling 
parameters flT5|) and flT6l) are nonvanishing, proceeds as follows. The method rests on the 
calculation of the mean-field energies of Bloch states Sm, which are obtained by numerically 
solving the GP equation (^ for the order parameters tprn |6|, in order to evaluate the integral 
yielding the energy per particle 

£m = J i^l'^^ral^ + ^^^p \^m\^ + ^N Q I^^M dx dy , (Al) 

where we may distinguish noninteracting and interacting contributions, 

C = / (l^l^^-l' + ^trap \^m?\ dxdy, (A2) 
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and 



C*= I^Ng\^^\'dxdy, (A3) 



respectively. In the context of this paper, i.e., for large barrier heights, the above GP 
energies ( ]A2l) -( !A3l) must coincide with those given by the corresponding contributions in 
f l?T]) with Jeff replaced from ( 1T^ . Then, by calculating (lA2p -( 1A3|) for any two Bloch states 
and equating such results to the corresponding terms in ( 12T1) . one can construct a linear 
system of four equations from which we may obtain the BH model parameters flT4|) - flT7|) . 
For instance, the simplest choice for Nc even corresponds to the winding numbers m = 
and m = Nc/2, which yields the following set of equations 

S^ = e- uJ, (A4) 

Slr'=^-^^{U/2-uJ'), (A5) 

£%./2 = e + uJ, (A6) 

^;7.% = ^^^(f//2 + ^J0, (A7) 

and the solution of such a system is given by the expressions fl2B]) - fl2^ . 
Finally, a similar calculation for N^ odd yields 

e = [^(V.-i)/2 + ^0 cos{7r/N,)]/[l + cos{n/N,)], (A8) 

2Nr. 



u 



= ]^[^S-i)/2 + ^rcos(7r/Ar,)]/[l + cos(7r/Ar,)], (A9) 

J = ^[^(V-i)/2 - ^o]/[l + cos(7r/Ar,)], (AlO) 

J' = 7T7#^[^(t^i)/2 - ^r]/[l + cos(7r/iVj]. (All) 



2{N - 1) [^(^-^/^ " ^° ]/f^ + cos(7r/iVj]. 
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